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SINGULAR VALUES FOR PRODUCTS OF COMPLEX GINIBRE 
MATRICES WITH A SOURCE: HARD EDGE LIMIT AND 
PHASE TRANSITION 

PETER J. FORRESTER AND DANG-ZHENG LIU 


Abstract. The singular values squared of the random matrix product Y = 
GrG-i —i---Gi(Go + A), where each Gj is a rectangular standard complex 
Gaussian matrix while A is non-random, are shown to be a determinantal point 
process with correlation kernel given by a double contour integral. When all 
but finitely many eigenvalues of A*A are equal to bN, the kernel is shown 
to admit a well-defined hard edge scaling, in which case a critical value is 
established and a phase transition phenomenon is observed. More specifically, 
the limiting kernel in the subcritical regime of 0 < 6 < 1 is independent of b, 
and is in fact the same as that known for the case 6 = 0 due to Kuijlaars and 
Zhang. The critical regime of 6 = 1 allows for a double scaling limit by choosing 
6 = (1 — r/y/N)~^, and for this the critical kernel and outlier phenomenon 
are established. In the simplest case r = 0, which is closely related to non¬ 
intersecting squared Bessel paths, a distribution corresponding to the finite 
shifted mean LUE is proven to be the scaling limit in the supercritical regime 
of 6 > 1 with two distinct scaling rates. Similar results also hold true for the 
random matrix product TrTr—i • ■ ■ Ti{Gq-\-A), with each Tj being a truncated 
unitary matrix. 


1. Introduction and main results 

1.1. Introduction. The squared singular values of a matrix X are equal to the 
eigenvalues of the positive semi-definite Hermitian matrix where X* denotes 

the Hermitian conjugate of X. An ensemble of random matrices of the form X*X 
may then contain x = 0 as the left boundary of support of the eigenvalues. Since 
the eigenvalue density is strictly zero for x < 0, x = 0 is then called a hard edge (see 
e.g. m Ch. 7]). As an explicit example, consider the ensemble of n x iV (n > N) 
rectangular standard complex Gaussian random matrices, namely the joint density 
of elements being proportional to exp{—tr(A*A)}, and let A be a matrix from 
this ensemble. Let {Aj} denote the eigenvalues of the scaled positive semi-definite 
matrix N~^X*X. In the limit TV —>■ oo with n — N fixed, the density of {Aj} 
has support [0,4]. That the support is a finite interval gives rise to this particular 
scaling being referred to as global scaling, and the corresponding density as the 
global density. The explicit functional form of the global density is given by the 
so-called Marchenko-Pastur law (see e.g. [54]) 

omW = ^ s 4. (Li) 
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Note in particular the reciprocal square root singularity as the hard edge A = 0 is 
approached from above, in contrast to the square root singularity as A —>■ 4“. The 
point A = 4 is an example of what is termed a soft edge, since for finite N the 
eigenvalue density is not strictly zero for A > 4. 

Continuing with this example, for large N the eigenvalues in the neighbourhood 
of the hard edge have spacing 0(1) upon the introduction of the scaled variables 
Xj = AN'^Xj (j = 1, N) (see e.g. [27l §7.2.1]). This will be referred to as hard 
edge scaling. Moreover, in the limit iV —>■ oo, and with vq = n — N, the limiting 
state — referred to as the hard edge state — is an example of a determinantal point 
process, meaning that the fc-point correlation function can be written in the form 

P(fc)(Xi,...,Xfe) = det[/C‘'(X„XO],,;=i,...,fe (1.2) 

with correlation kernel (see e.g. m Exercises 7.2 q.l]) 

1 

= J^oiV^)JuoiV¥)dt, (1.3) 

where is the Bessel function of the first kind of order z/q. 

Our interest in this paper is in the functional form and analytic properties of 
the correlation kernel for the hard edge scaling of the squared singular values of the 
product of independent random matrices 

y = G,G,_i---Gi(Go+Al), (1.4) 

where each Gj is an (iV + Vj) x {N + iyj-i) standard complex Gaussian matrix (also 
referred to as the complex Ginibre matrices since such non-Hermitian random ma¬ 
trices in the square case were first studied by Ginibre [34) 1 with u_i = 0 and integers 
uq, ..., Ur > 0, while A is of size {N + i/q) x N and fixed. If we focus on the singular 
values of Y, the definition CH) can equivalently be written as a product of inde¬ 
pendent square matrices Gj now with each being distributed according to the joint 
density of elements proportional to det*^^ (G*Gj) exp{—tr(G*Gj)}; see e.g. [51 H5]. 
In this case the restriction on the parameters can be relaxed to uq, ui,... > —1, 

and the main results in the present paper (for instance. Proposition 1 1.11 Theorems 
11.2111.3113. H and 13.21) also hold true. 

In the case that all entries of A are zero, the determinantal representation of the 
joint eigenvalue density and the limiting hard edge state have been the subject of 
a number of recent works [51 aiMi Ei m [57]. For r nonzero, the product (El 
is the simplest nontrivial example of the more general product {Gr + ^r) • • • (Gi -I- 
^i)(Go + 2 I 0 ), where each Gj is random and each Aj is fixed. However for the 
latter product with r > 1, it is not known how to find a closed form of the joint 
eigenvalue density, which is the starting point of our study (see Proposition 11.11 
below), let alone to study asymptotic statistical properties. 

The study of products of random matrices goes back to the pioneering work of 
Furstenberg and Kesten [33j in the context of dynamical systems and their Lya¬ 
punov exponents. Later, applications were found in Schrddinger operator theory 
m in statistical physics relating to disordered and chaotic dynamical systems [20] , 
in wireless communication networks [59] and in combinatorics [56] . Further motiva¬ 
tions from theoretical physics include the DMPK equation in mesoscopic quantum 
transport [TTJ|53|5T], fluid turbulence [551 eqn.(16)] and time evolution models [SOI 
eqn.(l)] (the last two models can be immediately recognized after discrete integra¬ 
tion in time). 
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Products of complex Ginibre matrices, and of truncations of Haar distributed 
random unitary matrices, have attracted much attention as examples of determina- 
tal point processes with kernels possessing special integrability properties. The first 
advance in this direction was by Akemann and coworkers IS [3], who derived the 
joint eigenvalue density and corresponding correlation kernel in terms of orthogonal 
functions for the product of complex Ginibre matrices; a double integral formula for 
the correlation kernel was subsequently obtained by Kuijlaars and Zhang, cf. |461 
Prop. 5.1]. These advances have opened up the possibility to study local statistical 
properties, see gsiiMiiisiin] for the hard edge limit and [ 3 S] for bulk and soft edge 
limits. All these studies form part of a fast paced and very recent literature relat¬ 
ing to the integrability and exactly solvable properties of random matrix products. 
Works relating to this theme which have appeared on the electronic preprint archive 
over the past few months (as of July 2015) include [48l [35l HOj IHl [39l | 6 TJ [32l [19]; 
we refer the reader to [5] for a recent survey article. Here we contribute to this line 
of research by undertaking a comprehensive study of the hard edge state formed 
by the singular values of (na. 

The first point to note is that the hard edge state in the case A = 0 depends 
on r, and thus is no longer described by the correlation kernel (11.31) . This fact can 
be anticipated by an analysis of the global density of the squared singular values 
[51 [5n[551[55] . The global density, which refers to the limiting density of eigenvalues 
of N~'"~^Y*Y as N ^ oo, is found to exhibit the hard edge singularity (see [56] or 
(an eqn (2.16)]) 

-sin^^A“^+^ as A ^ 0+, (1.5) 

TT r -|- 2 

which has an r-dependent exponent. In fact with the eigenvalues of Y*Y scaled 
according to Xj = Nxj (j = 1,..., Af), as iV —>• oo the hard edge state in the case 
A = 0 forms a determinantal point process with limiting correlation kernel 




p — 1 / 2 +ioo 


(27rf)^ 


- 1/2 — 200 


du (j) dt 




~Io 


0, -Uo, -Ui, . . . , -l/r 


r(uj -b u J- 1) sin?™ x*y 
r(uj -|- t -b 1 ) sin Tit u — t 


/ ’ V Uq , Ui, . . . , Ur, 0 


uy 


^ du, 

( 1 . 6 ) 


where denotes the Meijer G-function defined by the contour integral 

[ nr=ir(b,+g)n;^ir(i-a, -g) _ 

) 2TTi 4 rijLm+l r(l - bj - S) nj=n+l r(aj + sf 

(1.7) 

see [ini Sect. 5.2] for the choice of the contour 7 and elementary properties of 
G-functions, or [ 8 ] for a gentle introduction; it is worth mentioning a particular 
relation between the generalized hypergeometric function pFq and the Meijer G- 
function (cf. eqn (14), [49l Sect. 5.2]) 


a 


n,n ( 

V 


ai, 

bi, 


nLir(^z + l) ^i,p / l-ai,...,l-ap 
nLir(ai + l) P’^+H0,1-6 i,...,1-6, 


-z), ( 1 . 8 ) 
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which thus gives an alternative way of writing the first Meijer G-function in the final 
line of (jl.6l) . These kernels were described in [46] and are named after Meijer G- 
kernels in [45] . They also appear in the hard edge scaling for products with inverses 
of Ginibre matrices [29], products of truncated unitary matrices |45| . Cauchy two 
matrix models [T^ US] [30] , and Muttalib-Borodin biorthogonal ensembles [iniES] 
(cf. [45] for the relationship between Borodin’s expression and Meijer G-kernels). 

As noted in [46l Sect. 5.3], in the case r = 0 the facts that 


^ 1,0 

'-^ 0,2 


( 


0, —V 




{uyr/^M2^), ( 1 . 9 ) 


show 

y) = A{ylxr/^K'-{Ax, Ay). ( 1 . 10 ) 

The factor of {ylxyl"^ cancels out of the determinant (11.211 . while the factors of 4 
are accounted for by this same factor being present in the scaling leading to (11.31) : 
recall the text leading to this equation. 

Consider now m with 


A — ^/^I(N+uo)xN, ( 1 - 11 ) 

where I(N+vo)y.N denotes the {N A-vq) x N rectangular matrix with I’s on the main 
diagonal, and O’s elsewhere. It was shown recently in [an Remark 3.4] that there 
is a critical value of 5 = 1 for which as N ^ oo the left hand edge of the support 
of the global scaled squared singular values equals 0 for the last time as b increases 
from 0. Moreover, it was shown that the singularity of the global density has the 
leading form 

1 27r 2 

“ sin --_;v“i+ 2 P +3 as A —?► 0^, (1-12) 

TT 2r + 3 

which gives rise to a different family of exponents to those in (11.51) . We remark that 
the fractional part of the exponents, l/(r + 2) and l/(r + 3/2) respectively in (11.51) 
and ( 11 . 121 ) . are the reciprocals of positive integers and half-integers, which given 
knowledge of the correlation kernel (11.61) and its analogue in relation to ( 11 . 121 ) to be 
established herein (see eqn. (11.241) below), is coincident with them being the simplest 
in terms of tractability of the general rational fractional exponents accessible in the 
Raney family (see e.g. [ST] eqn. (2.16)]), so named due to the sequence formed by 
the moments of the global density. 

Let A be again given by (11.111) . and consider the case r = 0 in (11.41) so that 
T = Go -|- ■\/WV/( 7 v+i/o)xiv- It is well known that the squared singular values allow 
for an interpretation as the positions of non-intersecting particles on the half line 
evolving according to the squared Bessel process with parameter d = 2(uo -I- 1) 
(see e.g. EZllSH]). In this interpretation the particles all begin at the same point 
bN, evolve for t = 1 time units, and furthermore are conditioned to remain non¬ 
intersecting if the process was to continue to t > oo. The support of the density 
of such a process with a delta function initial condition NS{x — a) is, for Nt < a 
equal at leading order to [—+ a,Nt + a] (this fact is implied by results in [43], 
for example), so we see that with a = bN and t = 1, the particles first come in 
contact with the wall x = 0 as 6 is decreased to 6 = 1. A functional form of the 
hard edge scaled kernel in the critical case 6 = 1 , generalised to a double scaling by 
setting 6 = (1 — t/'/N)~^, has recently been obtained in [44]. In the present paper 
an alternative functional form to that in [44] is derived; see eqn. (11.251) below. The 
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kernel (ll.25|l . further specialised to uq = —1/2 reads 


1 

271^1 



g — TU+jU^+TV + jV^ 

u — V 


cos( 2 -yitJ) 


cos( 2 y^). 


(1.13) 


With ^ and rj replaced by squared variables, (11.131) is identified in [44] as the sym¬ 
metric Pearcey kernel found in the study [16] . Moreover, our method of derivation 
of this new functional form in the case r = 0 works equally as well for the double 
scaling of the critical kernel in the general r case, which is our main theme. The 
resulting explicit double contour integral expression is given in Theorem ll. 21 below. 


1.2. Main results. In preparation for the statement of our first key result, let us 
introduce two auxiliary functions. The first is defined to be 


4 '( m ; x) 


1 1 
(2717)’' r(i/o + 1) 


dwi 


dWr 'i 




1=1 


X (uo + 1; -ux/ (wi • • • Wr)), 


(1.14) 


where 71 ,..., y,. are paths starting and ending at negative infinity and encircling 
the origin once in the positive direction, or equivalently (cf. (|4.4|) . Sect. |4| below), 


1 If 

4 '( m ; x ) = J d,w{-x)~'^iFi{w;iyo + l;u) 

xr(R;)]^ ^ -r, (1-15) 

r(ui -I- 1 - w) 

with 7 encircling all non-positive integers, while the other reads 

-I /•c+ioo 

'^{v;y) = ^ ds 2 /"®(/( 7 ;;s)]^r(u/-f s), (1.16) 

J C —200 

where c > — min{uo, ui,..., Vj.}, and 

1 

= Fi —TTf / + (1.17) 

r(uo + 1) Jo 

^ r(!^o + 1 ) + s;vo + 1 ; -v). (1.18) 

In the case r = 0 (11.141) is to be interpreted as 

4'(M;a;) = oFi{i^q + 1;-ux), (1.19) 

and a calculation shows that (11.161) simplifies to read 

2 /) = y'''’e~y oFi{vo + 1; -vy). (1.20) 

The two auxiliary functions appear in a double contour integral expression for the 
correlation kernel, which we present next. Its significance is that it provides the 
starting point for further asymptotic analysis. The special case r = 0 was previously 
obtained by Desrosiers and one of the present authors; see [Ml Prop. 5]. 



















6 


PETER J. FORRESTER AND DANG-ZHENG LIU 


Proposition 1.1. Let Y be defined by and suppose that all eigenvalues 

ai,..., ON of A*A are positive. The joint density of eigenvalues for Y*Y can be 
written in the form 

Vn{xi,...,xn) = ■^det[KN{xi,Xj)]^j^.^ ( 1 . 21 ) 

with correlation kernel 

1 poo p -I ^ I 

KN{x,y) = — du dvu''°e~'^''~'"'it{u;x)^{v,y) - TT ^ , (1.22) 

Stti Jo Jc u-vfj^v + ai 

where C is a counterclockwise contour encircling —ai,..., —on but not u. 


Remark 1.1. When some of the parameters afs are null, the double integral repre¬ 
sentation ( 11 . 221 ) remains valid provided that du is interpreted as limg_,,o+ du, 
or for given it > 0 C is chosen such that Reju} < u with any v € C. 


One of the main results in the present paper concerns a double scaling limit near 
the critical point, which permits a new family of limiting kernels. 

Theorem 1.2 (Critical kernel). With the kernel (11.221) . for t G R let 

ai = ••• = OAT = 1V(1-T//iV)-^ (1.23) 


Then we have 


lim —=Kis[ 
N^OO y/Jf 


(_L JL) 




dv 




— TU—^U +TV+-^V 


U — V 


^ ^oJ-|- 2 ( 


0 ,-uo,-ui,.. 


— Ur 



^r+ 1,0 

^ 0 ,r +2 


(uo,Ui,...,Ur ,0 



=:/Ch'’-(^,ry;r), (1.24) 

valid uniformly for Cv 'in- any compact set of ( 0 ,oo) and for r in any compact set 
o/R. 

In the special case r = 0, upon making use of (11.91) we see from (11.241) that 


where the integral form on the RHS of the above equation is similar to ()1.22l) with 
r = 0 (cf. [Ml Prop. 5]). In the study [Ml displayed equation below (1.34)], this 
kernel was conjectured to be an equivalent functional form to that derived therein 
in the case r = 0. Our work thus provides a direct way of deriving ()1.25l) for the 
r = 0 critical kernel. Recently, an understanding of the resulting functional identity 
has been given in |22l Remark 2.26]. 

Theorem ll. 21 quantifies the limiting correlation kernel for the situation that Ok = 
N{1 — (fc = 1,..., N), which is shown to depend on r, thus justifying the 

term critical kernel. A variation on this setting is to have at most finitely many 
source eigenvalues, say oi,..., Om, go to infinity at a smaller but appropriate scale 
and others remain at the same critical value. This gives rise to a multi-parameter 
deformation of the critical kernel (|1.24l) . 


1 

27ri 


du / dv 


if) 


uo/2 ^-ru- 

U — V 

X (1-25) 
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Theorem 1.3 (Deformed critical kernel). With the kernel (11.2211 . for a fixed non¬ 
negative integer m, let 

Oj = '/Ngj, j = 1,..., m and Ok = N{1 — t/'/N)~^, fe = m + 1,..., IV, (1-26) 
where r G R and ai,..., > 0. Then we have 

^ Yj \ 1 r°° p—c+ioo 


u\^oe 

V y 


2m J-c-zoo 

m 

^^v + a-i ’ V0, —un, —ui...., —Ur 7 ’ V 


1 7^ + Ci 
1=1 J 


UO, Ul, ■ ■ ■ , Ur, 0 


.0, -Uo, -Ui,..., -Ui 

= : cr), 

where 0 < c < min{(Ti,..., am}- 

In the simplest case r = 0, upon making use of (11.91) we see from (11.271) that 


vqj 

(1.27) 





p—c+ioo 


27ri 


du 


dv 


uo/2^-ru-^u^+TV-i-^v'^ -^M + 


U — V 


n 

i=i 


V + (T, 


{2'\/vfi)J,jg (2y/vrj). 


(1.28) 


Even in this special case, the kernel (11.2711 appears to be new. 

We remark that the inter-relationship between the interpolating kernel (11.271) 
and critical kernel (11.241) is similar in form to that between the interpolating Airy 
kernel and Airy kernel (see e.g. [3 IB)- Furthermore, as the parameter b displayed in 
eqn. (11.1111 increases from zero, we will establish a phase transition at the hard edge 
from the Meijer G-kernel (cf. Theorem l3.1|l to the critical and deformed critical ker¬ 
nels (cf. Theorems II.21 and II.31) . then to the shifted mean LUE kernel (cf. Theorem 
13.21) : see Section [3] for more details. A similar phase transition occurs in another 
random matrix product T^Tr-i ■ ■ ■Ti{Go A), with each Tj being a truncated 
unitary matrix; see Section |4l 

The paper is organized as follows. Section [2] is devoted to the joint eigenvalue 
probability density function (PDF) and a double contour integral representation 
for the correlation kernel of the squared singular values of the product (11.411 . The 
proof of Proposition II.II will be given, and the formulas for the average of the ratio 
of characteristic polynomials and a single (inverse) characteristic polynomial are 
also derived. In Section [3] the hard edge limits of the kernel in different regimes are 
evaluated, which include the proofs of Theorems ll.2l and ll.3l Our methods are used 
to similarly analyse the product of r truncated unitary matrices and one shifted 
mean Ginibre matrix in Section |4l In Section [5] further discussions on asymptotics 
for large variables, and some open problems, are presented. 


2. Eigenvalue PDF and double integral for correlation kernel 

2.1. Correlation kernels. Consider (II.4|) in the case r = 0. Let xi,... ,xn and 
ai,..., Oat denote the eigenvalues of Y*Y and A*A respectively. It is well known 
(see e.g. [231 Prop. 5], [271 §11.6]) that the eigenvalue PDF of the random matrix 
Y*Y is an example of a biorthogonal ensemble m 

Qn{xi,...,xn) = ^det[77i(xj)]5=i det[^i(a:j)]5=i, 


( 2 . 1 ) 
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where r]i{x) = ^i(x) = x''° e~^oFiiv^ + l]aix), and Zn denotes the normali¬ 

sation. Our first task is to specify a functional form for the joint eigenvalue PDF of 
Y*Y in the case of general r. For this purpose use will be made of a recent result 
due to Kuijlaars and Stivigny |45| . 

Proposition 2.1 (Special case of [33 Thm. 2.1]). Let W be an n x n random 
matrix, and suppose that the eigenvalue PDF ofW*W can be written in the form 

{xk-Xj)Aei[fk-i{xj)]lf,^^ (2.2) 

l<j<fc<n 

for some {fk-i{x)}k=i,...,n- For v > 0, let G be an (n -\- v) x n standard complex 
Gaussian matrix. The squared singular values of GW, or equivalently the eigenval¬ 
ues of {GW)*GW, then have their PDF proportional to 

n (2/fc-yi)det[3fc-i(yj)]"_fc=i, (2.3) 

l<j<fc<n 

where 

9k{y) = [ x''e~'"fk(-)—, fc = 0,...,n-l. (2.4) 

Jo kx/ X 

Let Y be defined in (El and let ai, ... ,aN denote the eigenvalues of A*A. 
Starting with El, application of Proposition 12.H r times in succession shows that 
the joint eigenvalue PDF of Y*Y is equal to 

det[?7i(a;j)]5=i det[C*(a:j)]5=i, (2.5) 

where rii{x) = and (x) = $(—a^; x), while with T = ti - ■ -tr 

1 ro° j-oo r 

+ 1) y dti---J dtr'[[t'''~^e~*'{^Y°e~^oFi{i'o + l;-v^), 

(2.6) 

valid for r > 1 (for r = 0 ^j{x) is defined as below El)- Here ^{v,y) is actually 
the same as defined in (|1.16|) . for which application of the Mellin transform gives 

r 

y"~^<^{v, y) dy = (/)(-i;; s) F(u; -h s), (2.7) 

1=1 

while use of the inverse Mellin transform gives the sought expression. We stress that 
when some of the a^s in (12.511 coincide L’Hospital’s rule provides the appropriate 
eigenvalue density. 

The significance of the structure (12.51) is that it provides a systematic way to 
compute the corresponding fc-point correlation function P(fc)(xi,... ,Xk), where the 
normalization has been chosen such that integrating gives N{N — 1) • • • (TV — 
see e.g. m eqn (5.1)] for the definition. 

Proposition 2.2 ([151 Prop. 2.2]). With gij := r]i{x)^j{x) dx, let be 

invertible for each n = 1, 2,.... Defining aj by 

(k,]5=i)‘ = (bu]M=i)’', (2-8) 

and setting 

N 

KN{x,y) = Yi 

ui=i 



( 2 . 9 ) 
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we have that the k-point correlation function is given by 

P{k){xi,.--,Xk) = (2.10) 

We are now ready to complete the proof of Proposition 11.11 


Proof of Provo sition \l.l\. Our first task is to compute := Vi{x)fj (x) dx. For 
this purpose, we require the fact (see e.g. [U eqns. (6.2.15), (6.2.33) and (4.5.2)] ) 
that with L'f^(jj) denoting the Laguerre polynomial of degree n, one has the Hankel 
pair 


y poo 

L7{y)= / i""+”e-‘oFi(uo + l;-2/t)* (2.11) 

n!r(i/o + 1) Jo 

and 

n'e* r°° 

= FT-TTT / y""i«”(y)e-"oFi(po + l-,-ty)dy. (2.12) 

r(uo +1) Jo 

Combination of (12.71) , (11.171) and (]2.11l) shows that 

r 

5.,, = {t-iy.e^^L‘'f,{-a,)l[r{ni+i). (2.13) 

1^1 


According to Proposition 12.21 must now invert the matrix (|2.13p . With 
G = [gij]f^j^n let C = the entries Cij of G then satisfy 

N r 

e“'=^(*-l)!L('«i(-afc)nr(o+i)ci., (2.14) 

2=1 Z = 1 


Without loss of generality we assume that ai,..., gn are pairwise distinct. In this 
case the above equations imply 

Nr N 

^(*-i)!Lr^(«)nr(^'+*)"^^=e”“^ n ( 2 - 15 ) 

i=l 1=1 l=l,l¥=j 


as can be verified by noting that both sides are polynomials of degree iV — 1 in u 
which are equal at N different points since (I2.14|) is satished. Using this implicit 
formula for {cij} we now want to show that (12.91) implies the double contour integral 
formula (11.221) . 

Using the integral representation of the reciprocal Gamma function 


1 

W) 


A f w-^e^dw, 

2iii 


7 


(2.16) 
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we have from (j2.9ll that 


1 f f I I 

Kn{x, y) = ^j{y) / dwi--- 

i,,=l Jii Jir 




^ (zz~~~zr)' + 


"Wi ’ — Wr 
N 


Z=1 


( 0 ^-\r [ dwi--- f dwr'[[wi''‘ 

i=l "'71 Jlr 

poo 

duu''°L’^'!_i{u)e~'^oFi{i'Q +1; - — -) 

Jo Wi---Wr' 


r(uo + 1 ) 


1 f f 

dwi--- dwrl[wl 
J 71 hr 


-Vl-l Wi 


g glill-'-U 


poo 

/ duii'^“e-“o-Fi(uo + l;-e~°^ TT ~ ■ (2.17) 

Jo Wi---Wr aj - ai 


r(uo + 1) Jo 

Here the formulae (12.121) and (12.151) have been made use of respectively in the 
second and third equalities. 

Finally, with (12.61) and (I1.14|) in mind, these facts substituted into (12.171) imply 
that 


poo 

KN{x,y)= / duh°e~'^-^{u;x)Y,^i-(^j'^y) e"”" TT ~~ 

“'o i=i 

We recognise the sum over j as the sum of the residues at {a;} of 


ai 


(2.18) 


$(-?;;y)- 


1 


N 

n 


-u - ai 


V- ai 


(2.19) 


considered as a function of v. Applying the residue theorem and changing v to —v, 
we thus arrive at the desired result. □ 


Remark 2.1. The case in which each oj = 0 has been analysed previously [3l l4l 146) . 
but using different working. Thus instead of computing the inverse matrix (j2.8p . 
functions 


Pj-i{x) e Span{yi(x),... ,yj(x)}, Qj-i(x) G Span{^i(a;),... ,,fj(x)}, 

with the biorthogonality property Pk{x)Qi{x) dx = 6k,i were constructed. In 
terms of these functions (12.91) simplifies from a double sum to the single sum 

N 

KN{x,y) = Yh-i{x)Q,-i{y)- (2.20) 

1=1 


Instead of (11.221) with each ai = 0 (which strictly speaking is ill-defined due to 
the restriction on the contour C, but can be well understood in a limiting sense, 
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cf. Remark [ni) , this leads to the double integral formula 


KN{x,y) 


1 

(27ri)2 


^ — 1/2 + 200 


' — 1/2 — 200 


du ® dt 


i=-i 


r(j/j +1 + 1 ) 


r{t-N + l) 

T{u — + 1) u — t 


( 2 . 21 ) 


where E is a simple closed contour encircling anti-clockwise t = 0,— 1 
but not m; see [461 Prop. 5.1] for the detailed derivation, where similar integral 
representations for both multiple orthogonal functions Pj-i and Qj-i are first 
derived and then the double integral follows from a particular combination. 


Next, we further investigate Proposition 11.11 and establish a corollary under the 
assumption that all but a fixed number of source parameters are equal to a. Pre¬ 
cisely, for TO > 0 let ttm+i = ■■■ = ttN = a- More definitions are also needed. For 
= 1,2,3,... and n = 0,1, 2,..., set 





k-1 

duu''°e~'^'i>{u;x){u + a)" + o-i), 

1=1 


( 2 . 22 ) 


and 

k 

£^n\x;a,ai,...,ak) = [ dv e'"^{v;x){v + —, (2.23) 

2711 Jc zEi 

where C is a counterclockwise contour encircling —a, —oi,..., —Uk but not any point 
on the positive real axis. 


Corollary 2.3. Let Km be the kernel (11.221) . and for m > 0 let 

Hm-t-l = • • • = Om — u. (2.24) 


Then we have 


KN{x,y) = [ du [ dvu''°e “+''4'(ii; ®)$(-i;; y)^— 

2771 Jq 7c u — v\v + ay 

m 

+ X] a, ai,..., Ofc-i) cP_^{x; a, m,..., Ofc), 

k=l 

where C is a counterclockwise contour encircling —a but not u. 

Proof. We will use the identity 

1 ™ I 

1 TT 71-1-0; 


N-r. 


(2.25) 


n 




nf=i^(M + a;) 


(2.26) 


u-vf-J^v + ai ^ nti(^ + a;) 

which has been proved by induction in |23j : see the equation (5.12) therein. A 
direct proof can be given as follows. Rewriting u — v = u + ak — (v + a^), we have 


m -ttA;—1/ , \ m k 

(^_^)^ ^ ^ ,« + a; 


k-1 


fe=i n;’=i(^^ + «o 


fc=l 1^1 


V ai 


xn 

k=l 1=1 


u + ai 
V + ai 


n 


M -I- a; 
V + ai 


- 1 , 


(2.27) 


and (12.261) follows. 

Recalling (12.241) . substituting (12.261) in (11.221) . comparing the sought equation 
with (12.221) and (12.231) this completes the proof. □ 
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2.2. Average of characteristic polynomials. Recall that a biorthogonal ensem¬ 
ble [15] refers to the joint density function 

Qn{xi,...,xn) = ^det[77i(xj)]5=idet[Ci(a;j)]5=i, (2.28) 

where all variables Xi,... ,xn are assumed to lie in the same interval ICR for 
simplicity. 

For the special case rji = the average ratio of characteristic polynomials 

under the density (I2.28P can be expressed in terms of the correlation kernel; thus 
as a minor variant of |24l Prop. 1] we have the following. 


Proposition 2.4. With the same assumption and notation as in Provosition \2.Si 
let r)j{x) = x^~^ for j = 1, 2,.... Then, for z € C\I 


E[TT^^] = [ du- — -Kn{x,u). 


Equivalently, if for a; € R ire define the residue 


Res 2 =a;/(z) = lim -lmf{x - ie), 

£-> 0 + TT 


then 


K]sr{x,y) = ^^Res^=a;E[ 


N 


x-y 


1=1 


X- Xl - 

z- xr 


(2.29) 


(2.30) 


(2.31) 


In the case of the average of a single characteristic polynomial or its reciprocal, 
alternative expressions are also available; cf. Proposition 2 of 


Proposition 2.5. With the same assumption and notation as in Provosition \2.2\ 
let r]j{x) = x^~^ for j = 1,2,.... Then, 


N 


^fj^z-xi^ Zn 


91,1 


91,2 


9n-i,i 9n-i,2 
f du ^ Ldu^ 

JI Z — U JI z — 'i 


f 1 
Jj 


for z € C\I and 


N 


m 


E[J|(a:-Xi)] = ^ 


z=i 


91.1 

92.1 



9i,N 


9n-i,n 

) 

i 


91 ,N 

Vlix) 

92,N 

V2ix) 

9n+i,n 

VN+l{x) 


A-l-l 


CA+1,A-|-1 




N+lX‘ 


J-l 


N 


(2.32) 


(2.33) 


(2.34) 


(2.35) 


where the normalization Zjq = N\ and Cj^N+i is the {N + l,j) entry 


of the inverse of 
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Proof. The formulas ()2.32|1 and ()2.34|1 have been proved in Proposition 2 of |24j . 
After noting the facts [gi,j]N[cij]% = In and N\/Zn = det[cij]^, building on (12.3211 
simple manipulation gives (12.331) . On the other hand, by [ci^jY'Nj^i[gij]N+i ~ ^n+i 
and 

N\ (iV+l)! Zat+i u 

Zn {N + 1)Zn Zn+1 {N + 1)Zn^ (• ) 

we have from (|2.34p that 


N 






iV+l 


N+1 


{N + 1)Zn 




N+lX^ 


(2.37) 


1=1 ' j=i 

which further implies the sought equation p.35|) since it is a monic polynomial. □ 

Application of the previous two propositions gives us explicit evaluation of aver¬ 
ages of characteristic polynomials for the product of random matrices m- 

Proposition 2.6. For the eigenvalue PDF i2.5\} . the following hold true. 

(i) Let Kn be the kernel (11.221) . then for z £ C\IR, 


N 


E[n 




= / du^KNM. 
z-xp Jq z-u 


(2.38) 


(a) Let $ be given by (11.161) . then for z G C\R, 


N 


E[n 


(- 1 ) 


N-1 


1=1 


z- xY T{ui + N) Jq 


pOO p 1^1 


where C is a counterclockwise contour encircling —oi,..., —on but not any point 
on the positive real axis. 

(Hi) Let di be given by (11.141) . then 

N r ^ N 

+ / duu''°e~'^'i!{u]x)'^{u + ai). (2.40) 

i=i 1=1 1=1 

Proof. It is immediate that Proposition 12.41 implies (i). For (ii), noting that the 
leading term of the Laguerre polynomial is 


n\L'0[x) = {-xY + ■■■ , 

dividing by and taking the limit u —>■ oo in (I2.15|l we see that 

N 




n 


1 


n 


i^Ar(uz + iV) 


(2.41) 


(2.42) 


Substituting cnj in (I2.33|l and noting gj{u) = $(—aj;u), we obtain (12.3911 . 

For (iii), we first introduce an auxiliary variable on+i and set rjN+iiu) = 
^{-ON+pu). The fact that {[gij]N+iy[cij]N+i = In+i implies 


cat-i-iw-i-i = 


det[g^j]Af 

det[gij]Ar+i r(u; -f A^ -f 1) 


n 


-O^N + l 


N 

n 

1=1 


1 


Qn+i — ai 


(2.43) 
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Changing TV to + 1 and using (I2.15I1 . as derived in (j2.17|) we obtain 

N-\-l N 

Cj TT — - —(2.44) 

^ ’ Jo o.N+1 - ai 

Combination of (12.431) , (I2.44p and (12.351) completes the proof of (iii). □ 


Again, for the eigenvalue PDF (12.5|) . let 


N 


Qn-i{x) = Res2=2,E[]^ 


1=1 


z - xr 


X £ (0, oo), 


then use of Proposition 12.61 (ii) shows 

1 (- 1)^-1 


N 


QN-lix) = 


2m ULin^i+N)Jc 


dv e’'4>(u; x) 


1 


V + ai' 


(2.45) 


(2.46) 


when oi,..., Oat are pairwise distinct it is a special case of Proposition 2 [M]. Also, 
let 

N 

Pjv(x) =E[[](x-t,)], (2.47) 

1=1 

combining Corollary 12.31 where m is taken to be zero and Proposition 12.61 the 
correlation kernel ATat given by (11.221) can be expressed as the single sum ()2.20l) 
in terms of Pj{x) and Qj{x). Here, without loss of generality, it is assumed that 
Pj{x) corresponds to the multi-parameters ai,...,aj while Qj(x) corresponds to 

Oi, . . . , Oj+i. 


Remark 2.2. In the special case r = 0, use of (I1.19|) shows that (I2.40p reduces to 


Pn{x) 


(-l)^e^ 
r(uo -k 1 ) 


poo 

Jo 


N 

oFi{i^o + 1; -xu)Y[{u + ai) du. 


This same expression has been derived using combinatorial means in [28] , and as the 
solution of a partial differential equation in M- Furthermore, in this case Qn-i{x) 
and Pn (a;) are so-called multiple functions of type I and II respectively, and (12.201) 
reduces to Corollary 7 in [24]; see [24] or [42] for more details, especially when the 
parameters a^’s coalesce into D different values. For the case of ai = • • • = on = 0 
and general r, Qn-i{x) and PNix) are also multiple functions of type I and II 
associated with r + 1 weights; see [JS]- However, in the general case it remains as 
a challenge to identify a multiple orthogonal functions structure. 


3. Hard edge limits 

In this section we choose the source A such that all but possibly a fixed number 
m of the eigenvalues of A* A are equal to 6TV. Three regimes are distinguished: 
subcritical regime 0 < 6 < 1, critical regime 6=1 and supercritical regime 6 > 1; 
as to the former two regimes, see [23] and [22] for a relevant discussion on non¬ 
intersecting Bessel paths which corresponds to the case r = 0. In the present paper 
we focus on the scaled hard edge limits in the three regimes and leave the bulk 
and soft-edge limits to a future work; for the case oi = • • • = ajv = 0 the latter 
two limits have been established in [21]. The critical kernel results from a double 
scaling limit, and its functional form is our main result as stated in Section [l] As 
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b increases from zero, we will describe a phase transition from the Meijer G-kernel 
CH) to the critical kernel (cf. Theorem If .211 . then to the shifted mean LUE kernel 
(11.221) (cf. Theorem 13.21 for the case r = 0). 


3.1. Limiting kernels. We first suppose that 0 < 6 < 1. The hard edge scaling 
in this parameter range is in fact independent of b, and the hard edge correlation 
kernel (11.61) already known for the case 6 = 0 is reclaimed. 

Theorem 3.1 (Subcritical regime). With the kernel (ll.22[) . let 

oi = • • • = Oat = bN. (3-1) 

Then for 0 < b < 1, we have 

{l-b)N^^({l-b)N' (l-6)fv) 

where is given by ms), valid uniformly for r] in any compact set of (0, oo). 

Proof. Introducing rescaled variables x = ^/((l — b)N)^ y = rj/{(\ — b)N) in (jl.22l) 
and substituting u, v by uiV, vN respectively, we obtain 


(1 - b)N 


Kn{^ 




V 


{l-b)N' {l-b)N 




2(1 — b)TTi 


{NuY°'^{Nu; 


dv 


(1 - b)N 




(1 - b)N 


), (3.3) 


where 

f{z)=z-log{b + z). (3.4) 

Although both the functions ik and $ in the integrand of (13.31) depend on N, we 
will see that for the large N they do not enter the saddle point equation (e.g., cf. 
(I3.10|) and ()3.11|1 below). So we may perform saddle-point approximations and this 
is what we will do next. 

Consider now the exponent on the RHS of (13.31) . Since 


/W = l- 


1 

b + z' 


(3.5) 


there is a saddle point zq = 1 — b. We hereby deform the contour C into the union 
of two closed contours Ci {jCf such that Ci = {z & <C ■. \z + b\ = 1} and Cf is a 
clockwise contour encircling the segment [0,1 — 6 ] but not —6 {Cf and C 2 refer to the 
same curve except that the former indicates the clockwise direction). For instance, 
we can choose C 2 as the union of two segments from —0.56 to — 6 -l-e^*'^ respectively 
and an arc {z : z = — 6 -|-e®®,—e < 9 < e} for some small positive e. With such 
a choice, we divide the integration over C into two parts, and furthermore rewrite 
the double integral on the RHS of (13.3p as a sum of two integrals, that is. 






(l-6)iV "V(l-6)iV’(1-6)1V 


1 = p.v. / du di;(-)-I-p.v. / du dv{-) 
' Jo JCi Jo Jcz 

(3.6) 


:=/i 


10 

l2- 


Here the notation p.v. denotes the Cauchy principal value integral. It is worth 
stressing that, from (13.3L we can put some restrictions on the range of u,v in the 
above integrals such that u ^ 1 — b and v ^ ±i. This is done for the convenience 
of subsequent asymptotic analysis only. 
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As iV —^ c», we claim that the leading contribution of the double integral on the 
RHS of (j3.3|l comes from the range of u G (0,1 — 5) and u G C 2 . Actually, for / 2 , 
when u > 1 — b the u-integral vanishes by Cauchy’s theorem since the integrand 
does not have any singularities inside C 2 , while for 0 < u < 1 application of the 
residue theorem gives 

rl-b 


= 


1 


1-6 


du {NuY°'if (^Nu; 




(1 - b)N 


)i>{Nu- 


V 


(1 - b)N 


)■ (3.7) 


Using the asymptotic expansion of the function iFi for the large argument 
(cf. Theorem 4.2.2 and Corollary 4.2.3, [7]), for large N we have 


ifi(uo + s; uq + 1; -Nv) = {Nv) 

r(l - s) 


(l + 0(^)), Rei;>0, (3.8) 


and 

ifi(uo + 5 ; uo + 1; -Nv) = + ^(l)), Reu < 0. (3.9) 

r(j/o + s) A' 

Keeping in mind (11.1611 and (jl.l8|) . by definition of the Meijer G-function p.7D we 
have from (j3.8p that 


{Nu)''°i>{Nu; 


V \ ^r+l,of 


Vo, 


■ •, Vj -, 0 


ur] 


1-6 


(3.10) 


Here and below we use the notation /at ~ gjsr to mean that limAr_>.oo fN/gN = 1- 
On the other hand, consideration of the definition (I1.14|) shows 


'if{Nu;' 




-)-Go°+2{ 


0,-iyo,...,-Vr 


(1 - b)N 

where use has been made of the identity (cf. eqn (11.81) ') 

1 




1-6 


). 


(3.11) 


^ 1.0 / 

'-^ 0 ,r+ 2 \^ 


0, -Z/Q, ■ • ■ ,-Vr 


n 


r(ui +1) 


Fr+i(vo + 1 ,... ,1/^ + 1]—zy (3.12) 


Combining (13.7L (13.101) and (I3.11L and changing variables we get 


(3.13) 


Next, we deal with the integral Ii and show that it is negligible. In this case 
because of different asymptotic forms of iFi given in (|3.8D and (13.9p . we divide Ii 
into two parts as 


h = 


poo n poo p 

p.v. / du dv{-) + / du dv{-) := III + 112 , (3.14) 

Jo Jci,+ Jo Jci,-^ 

where Ci,+ = Ci(^{v : Reu > 0} and Ci,_ = Ci{^{v : Reu < 0}. Notice that 
for 0 < 6 < 1 one can easily check that Re{/(M)} attains its global minimum at 
u = 1 — 6 over (0, 00 ), while Re{/(u)} attains its global maximum at u = 1 — 6 over 
Cl. Therefore, for In combining (11.161) . (ll.lSp . (13.81) and (13.lip we have 


hi ^ p.v. 


du 


dv ^ e-^(7U)-/C)) u 

2(1 — b)Tri u — V 


G, 


1.0 ( 

0 .r+ 2 \^ 


0,-Vo,...,-iyr 


U^ 


1-6 


^r+ 1,0 / 
^ 0 ,r +2 


Vo,...,Vr,0 


VTj 


1-b 


(3.15) 
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For this, the standard steepest descent argument shows that the main contribution 
comes from the neighbourhood of the saddle point zq = 1 — b, namely, 

Ill = 0{1/Vn). (3.16) 

Similarly, for I 12 combination of (I1.16|) . (I1.18|) . (13.91) and (13.111) then gives us 


I 12 ^ du dv— -- . 

Jo Jci _ 2(1 - bjni u-v 


,^1,0 

Lzn 


^(, 




i-b 


X ■u'^Og-^f(log(b+«) + l-b). 


2717 


nc+ioo r 

^ Jc—ioo i 0 


) 




(3.17) 


for which the integrals of u and v respectively afford us bounds 0{l/^/N) and 
(^(^i.o+ 2 c-ig-(i-b)Ar)^ Together, we obtain the exponential decay estimation 

/12 = 0(7V'^o+2c-3/2g-(l-b)Ar)^ 

Combining (I3.13p . p.l6l) and (I3.18p . we arrive at the equation (13.21) . Further¬ 
more, it is clear that the previously derived estimates are valid uniformly for r] 
in any given compact set of (0, 00 ). □ 


Remark 3.1. When all the parameters aj’s are null, if we understand the double 
integral representation (11.221) as described in Remark ll.il then the same argument 
as in the proof of Theorem 13.11 is also applicable. This gives another derivation of 
(11.61) different from that in [46] . 

We turn to proofs of Theorems 11.21 and 11.31 


Proof of Theorem \1.2i Rescaling variables in (11.221) . we have 


du dv - 

2717 Jo Jc 77 - 7 ; y/N y/N 


(3.19) 


where f{z) = z — log(l + z/b) with 6 = (1 — r/y/N) If b is equal to the critical 
value 1, then the saddle point of f{z) is zg = 0. This time, for a fixed small number 
(5 > 0 we choose the contour as 


C = {z = -! + (! + 2(5)e*® : ffg < l^l < tt} U £a_oa+, (3.20) 

where 

do = arccos ^ ^ , A± =-1 + (1 + 2S)e^^^'’, (3-21) 

1 -b 2o 

and £a_oa+ denotes the union of two line segments from the point A_ to the origin 
to the point A_|_. It is clear that A± = (d, ±y /(2 + 3(5)(5), and the intersections of 
the 7/-axis and the contour C are B± = (0, ±2-\/(l + d)d). Moreover, the four 
points come close to the origin as d —> 0, which permits us to use the Taylor series 
expansion of f{v) for any v G C+ defined below (13.221) . 

First, we divide the integral on the RHS of (I3.19P into two parts 

1 (A _ 

vn 


du 


dv{-) 


du 


du (•):=/_+J+, (3.22) 


JC- 
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where C_ = {v € C : Rei; < 0} and C+ = {w € C : Rei; > 0}. We claim that 
the dominant contribution to ()3.19|1 comes from the neighbourhoods of ito = 0 and 
Vo = 0, so we need to expand the function f{z) at zo = 0. With the double scaling 
in mind, we obtain the Taylor series 


/W 


TZ 1 , T,,n 1, T n 

—F= + 1^(1- 7=) ^ o (1- 7=)^^^ 


(3.23) 


Therefore for /+, combining (I1.14|) . (I1.16|) . (11.181) and (13.81) . together with the 
relation (I3.12|) and the definition of Meijer G-function (jl.7|) we see that 


y/N 

2171 


du 


'0 


p-N{f{u)-f{v)) 

dv- --r“ 

U — V V 


^1,0 / 
^0,r+2 I 




^<)gZ+2{ 


UO, . . . ,Ur,0 


y/Nvri^. (3.24) 


Fix the two endpoints B± of C+ and deform it to the imaginary axis, then after 
substituting (I3.23P into (I3.24|) and rescaling u, v by u/y/N, v/y/N, we conclude that 
/+ converges to the kernel defined by (ll.24p . uniformly for ^,r] in a compact set of 
(0, oo) and for r in a compact set of R. 

Secondly, for the integral combination of (I1.14L (11.161) . (11.181) and (13.91) 
yields 


/_ 


27ri 


du dv- 


o-Nf{u)-N log(l+'!)/b) 




1.0 ( 


y/NuZ 

0, -Uo, ■ • ■ , I 

-I ^C+200 ^ 

— dsr]-%-vy-^N''°+^^^-^'>GY[r{izi+s). (3.25) 

ZTTI j c—ioo 7_T 


Since for sufficiently large N, 

Re{log(l + ^)} = i log ((1 + 26y + > log(l + S) (3.26) 

holds true uniformly for r in a compact set of R and for w G C_, use of the steepest 
descent argument leads to an exponential decay 


/_ _ (^^jyi'o-l+1.5Cg-Alog(l+(5)^ ^ 


(3.27) 


Lastly, by combining the foregoing results for /_ and /+, we then complete the 
proof. □ 


Proof of Theorem \1.3\ Rescaling variables in (11.221) . we have 


1 

y/N 



J. _ V_\ 

y/N' y/n) 


1 

27ri 



dv iV'^o + l/ 2 yivo 


,-N(f{u)-f{v)) 


U — V 


X 


n 


i=i 


{v + b){u + aj/y/N) 
{u + b){v + aj/y/N) 




(3.28) 


where f{z) = z — log(l + z/b) with 5 = (1 — T/y/N)~^. 

Proceeding as in the proof of Theorem II.21 Taylor expanding f{z) at z = 0, and 
rescaling m, v by u/y/N, v/y/N, we can complete the proof. □ 
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We next consider the supercritical case, that is 6 > 1. For r = 0, the limiting 
eigenvalue density has support [Li, L 2 ] with Li > 0 (thus the left-most end changes 
from the hard to the soft edge as b increases beyond unity as already remarked in 
the Introduction); see e.g. [43l|44]. However, for r > 0 and fixed uq, , Ur > 0, 
considerations from free probability theory suggest that the support will include 
the origin for general b. Nonetheless, in the simplest case of r = 0 a particular 
tuning and scaling of the supercritical case can be given which, on an appropriate 
length scale, effectively separates a bunch of eigenvalues near the origin from the 
rescaled left-end support. A similar result is conjectured to be true for the general 
r > 0. 


Theorem 3.2 (Supercritical regime for r = 0). With the kernel (jl.22ll where r = 0, 
for a fixed positive integer m let 

Qj = Ojb /(6 — 1), j = 1,..., m and Ok = bN, k = m + 1,..., N, (3.29) 
where 6 > 1 and cti, ..., am > 0. Then we have 


< / du I dvoTili^o + l;-uOoFi(i^o + 

Jo Jc 


1 


r] 


2TTi (r(z/o + 1))^ 


^—u+v 


n 


U + (T,- 


u — V V + aj 

j=i J 

where C is a counterclockwise contour encircling —ai ,..., —am but not u. 
Proof. Set K = [b — l)/b. For the large TV, we have from (11.221) with r = 0 that 


(3.30) 


(i-i)K„((l-i)Mi-i),) 


K 

2711 



dv u'^ 


^—u+v 


U — V 


X '^(u; K^)$(u; KT]) 

i=i 


u + aj/n 

V -I- aj!K 


/l + u/{bN) 
\l + v/{bN) 


N—m 

:=/i+/2, 


(3.31) 


where we have rewritten C = C 1 UC 2 and Ij = du dv (■). The closed contour 
Cl encircles and on its left lies the path C 2 encircling — &iV, 

beginning at and returning to — 00 . Keeping (ll.lOp and (11.201) in mind, for the 
choice of the contour of I 2 we have used the asymptotic property of qF’i (cf. [J^l 
Sect. 5.11.2]) 


1 N r(Uo + I) , 
oFi{vo + i;^) = „ ^ (-z) " 

Z^/ TT 

X (e-2*'^(l + 0(^))+e2*^(l + 0(^))), Izl^oo, (3.32) 

Again with (I1.19I) and (ll.20|) in mind, by taking the limit and changing variables, 
as —>■ 00 it is less difficult to know that goes to the desired integral. 

For the part I 2 , by the fact (I3.32I) , taking the limit in the integrand we see that the 
u-integral over the closed contour C 2 equals zero since the resulting integrand has 
no pole. The proof is thus completed. □ 


By comparison with (I1.22I) in the case r = 0, N = m, {ai} = {ai}, substitution 
of (I1.19I) and (ll.20p shows that the RHS of (|3.30l) is equal to e^~^Km{£,, ^)l{ai}={cri}, 
which is equivalent to the kernel for the m x m Laguerre Unitary Ensemble with a 
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source (see e.g. [24] or m Chapter 11]). For general r > 1, as to the supercritical 
case of 6 > 1 the following similar result is expected to be true 

N^oo^ b’ V b’’) 2Tri Jq Jc u-v 

™ -I- ■ ~ 

X 4'(u/it;K^)$(u/K;K7?) =: o'), (3-33) 

fJl ^ 

where k = {b — l)/b, di, $ are given by (11.141) . (I1.16|) . and C is a counterclockwise 
contour encircling —cri,..., —am but not u. To prove it, if we might control the 
behavior of k^) and e'^’'<l>(u; kt]) (for instance, we can try to derive an esti¬ 
mate KT]) = 0{v~^) as u —>■ —oo which should further be expected to vanish 

sub-exponentially), then /2 —>■ 0 and Ji goes to the desired integral as in the proof of 
Theorem 13.21 But such an estimate is yet to be found. Furthermore, we expect the 
correlations implied by (13.331) to be the same as those for Km{^,'ri)\{ai}={(Ti} after 
being multiplied by the factor g{K',r])/g{K',^) for some properly chosen function g. 
However, the mechanism which makes this true in the cases r > 1 remains to be 
clarified. 

Remark 3.2. It is of interest to contrast the scalings of {aj} in Theorem 11.31 and 
13.21 applying to the critical and supercritical cases respectively. Some insight as to 
the chosen values is possible by restricting attention to the case r = 0, for which 
the squared singular values have the interpretation as non-intersecting Brownian 
particles confined to a half line, as mentioned in the Introduction. In this interpre¬ 
tation, the initial position of particle j is aj, and the particles evolve for time t = 1. 
We interpret the values in aj in Theorem 11.31 as being such that the particles at 
the hard edge are all of the same order, with the k outlier particles appropriately 
merging with the spectrum edge of the N — k particles which started originally at 
iV(l — T/y/N). On the other hand, in Theorem 13.21 only the k particles starting at 
order unity from the hard edge are at order unity from the hard edge when t = 1, 
with the remaining N — k particles never reaching the hard edge by this time. 

Remark 3.3. If we strengthen the results in Theorems 11.2111.3113.11 and 13.21 from 
uniform convergence into the trace norm convergence of the integral operators with 
respect to the correlation kernels, then as a direct consequence we have the lim¬ 
iting gap probabilities after rescaling, especially including the smallest eigenvalue 
distribution; see Ea Chapters 8 & 9]. Since the proof of trace norm convergence 
is only a technical elaboration that confirms a well-expected result, we do not give 
the details. 


3.2. Characteristic polynomials. In this subsection we want to evaluate scaling 
limits for the ratio of characteristic polynomials according to three different regimes. 

Theorem 3.3. With the eigenvalue PDF i2.5\) . fix m € {0,1,2,...,} and let 


a-m+i = • • • = Cat = Nb. 

(i) Set Qj = Nbj with bj > 0 for j = 1,..., m, if 0 < b < 1, then for ( € 

N 


1 


.NAT 

N—^oo (1 — b)N 


E 


[n 


xi-f/i{l-b)N) 


xi-C/i{l-b)N)i 


du 




C-u 


(3.34) 
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(a) Set Qj = \/N(7j with aj > 0 for j = 1,... ,m, if b = 1/(1 — r/yfN) with 
r £ R, then for f £ C\R, 


lim — 



Xi - f/VN - 

xi - C/y/N- 



duj — 

C-u 


(3.35) 


(Hi) Set Gj = cFjb/{b — 1) with aj > 0 for j = 1,... ,m, ifb>l and m>l, then 
for r = 0 and for ( £ C\]R, 


lim 

N—^oo 


1 ^ 

(i-j)E[n 

1^1 


xi - (1 
Xi - (1 




du^}C)fi^,u;l 
c, u 



(3.36) 


Proof. By Proposition 12.61 Theorems 13.11 11.31 and 13.21 imply the sought results 
although a minor modification in the proof of Theorem l3.11 is required in relation to 
(13.341) (in the same circumstance the limiting subcritical kernel still holds true). □ 


Likewise, based on Proposition 12.61 we can prove the following theorem con¬ 
cerning the average of one single characteristic polynomial or its inverse. For this 
purpose we introduce four sets of generalised multiple functions (we say generalised 
since only for r = 0 do we know the multiple polynomial system; recall Remark 
12.21) of types II and I with m parameters cti, ..., am > 0. For k = 1,2 ... ,m, we 
define two sets of generalised multiple functions by 


poo 

FW(i;ai,...,afe_i) = / 

Jo 


X G. 


1,0 






k-l 


XU 


n(^^ + ^,), (3.37) 

i=i 


and 


^ -I pioo 

r^'^\x;ai,...,ak) = — / dvv-’'° 


^r+1,0 / 
^ ^0,r+2l 


vo,iyi,...,Vr,0 


IT-, (3.38) 

j=l J 

while for 0 < K < 1 two sets of Laguerre-like generalised multiple functions are 
defined by 

/■OO 1 

C^^\x-,K;ai,...,ak-i)= / du{u/e~^'l!{u/ k] kx)Y\{u + ai), (3.39) 


and 

fc 1 

(3.40) 

7 7 —I— m 

7=1 

Here 7 is a closed path which is encircling —a\,..., —Om once in the positive di¬ 
rection. 


k 

L^^\x-,K-,ai,...,ak) = J— [ dve'"^{v/K-,Kx)T \— 
2tti j , , ^ + 


Theorem 3.4. With the eigenvalue PDF \2.5\) . fix m £ {0,1,2,...,} and let 

Gm+l = ■ ■ ■ = ON = Nb. 
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(i) Set Qj = Nbj with bj > 0 for j = 1,... ,m, if 0 < b < 1, then 

N 


-Vn 

^Too ^5^ Lli XI - C/((l - b)N) 


[n: 


du 


1 ^r+l,o[ 


-G 


(-U °’"'+^Vuo,Ul,---,Ur,0 


i), (3.41) 


lim VNY[{iyi+N) 

1^1 1^1 

■i 


_ r^i.o 


f) (3.42) 


where 

r m 

T^ub) ^ (_i)Af^jV'^o+Afg-(i-b)Af {l-byo ]^(1 -b + bj). (3.43) 

1=1 1=1 

(a) Set Gj = '/Nuj with aj > 0 for j = 1,... ,m, if b = 1/(1 — r/yfN) , then 

N 


[n 


1 


-Vn^ 

lim ——r-E I , , , _ 

Llla;i-C/ViV 


r du^ 

lo C-u 


rW(M;i7i,...,a™), (3.44) 


r N 

=r(’"+')(e;ai,...,a„) (3.45) 

Z=1 /=1 


N—¥oo 

where 


= (-1)^N^+Go-ra)/2^y/Nr+T^Jr N). 


(3.46) 


Z=1 


lim 


(in) Set aj = ajb/ {b — 1) with aj > 0 for j = . ^m, if b > 1, then for r = 0, 

N 

E[n- 

“-t V- 0 ^ \J 3 

(3.47) 


-b 


1 


(6 - 1)T^"P) Ll^i XI-{1- i)CJ 


/ 


du- - 

0 C-u 


1-1/6; CTi,.. 


N 


ITGi (xz-(l-l/6)C) = /:(™+^^(C; 1-1/6; cti, ..., cr^) 

(3.48) 


/ = 1 


/ = 1 


where 


r>(sup) 

- N 


= (-l)^(6iV)"-™(6/(6- l))”". 


(3.49) 


Proof. By Proposition l2.61 following almost the same procedure as that in Theorems 
I3.11ll.3l and l3.2l we can evaluate the scaling limits. As a matter of fact, the proof will 
be simpler since it only involves a single variable integral. We omit the details. □ 


Let us conclude this section with two relationships between the limiting kernels 
(cf. (I1.27() and (j3.33p l and the generalised multiple functions defined by (I3.37p - 
(I3.40p : cf. Corollary 12.31 
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Proposition 3.5. We have 

m 

y; T, a) = y- t)+ {x; cri,..., cife-i) (y; cti, ..., Cfe), (3.50) 


y; ct) = ^ (x; k; o-i,..., ak-i) k; cti, ..., cTfc). (3.51) 

fc=i 

Proof. By use of the relation (I2.26L noting the definition of involved functions 
(13.37|) (13.401) . term-by-term integration immediately implies the above two formu¬ 
las. Here use has been made of /Cq’^(x, y; k; cr) = 0 for the second formula. □ 


4. Product with truncated unitary matrices 

The derivation of the double contour integral expression (11.2211 for the correlation 
kernel is expected to be applicable to a wider class of biorthogonal ensembles, 
specifically to those characterized by the form of (EU with rii{x) = X® ^ and £,i{x) = 
h{ai, x) for some appropriate function of two variables h and N generic parameters 
oi,..., ajv. In this section we consider the specific case of the biorthogonal ensemble 
corresponding to the product of r truncated unitary matrices and one shifted mean 
Ginibre matrix and derive a double integral representation of the correlation kernel 
and analyze the scaled limits at the hard edge. Other types of products ■ ■ ■ XiZ, 
where each Xj is a Ginibre or truncated unitary matrix while Z is a spiked Wishart 
matrix of the form GqS or a triangular random matrix (cf. |181132j b are presently 
under consideration [47]. 

Explicitly, instead of (lEl, we now consider the matrix product 

Y = Tr---mGo + A), (4.1) 

where each Tj is an (N + r'j) x (N -|- u^-i) truncation of a Haar distributed unitary 
matrix of size Mj x Mj and Gq is an {N + vq) x N standard complex Gaussian 
matrix while A is of size {N + i/q) x N and fixed. Here u_i = 0, vq, ... ,Vr are 
the nonnegative integers and y-j := Mj — N > Vj (for the general Vj > —1 the 
analysis below is also applicable). In the case that the matrix (Go -I- A) is absent, 
this product has been studied in a recent paper [40] . 

An analogue of Proposition 11.11 for the correlation kernel can be given. As 
in Proposition 11.11 two auxiliary functions are required, and so as to stress the 
structural similarities, analogous notation is used. Specifically, with r = 1,2,..., 
and 0 < g < r, the first is defined to be 


x) = 


(27ri)’' r(uo -b 1) y(o,oo)<* Jt 


dt / 




1^1 


X n exp { X ( 




Wi - • - Wr ■ 


1^0 + 1;-XU— - —Y (4.2) 

Wl-■ -Wr' 


where P = 71 x • • • x 7,., and 71 ,..., 7^ are paths starting and ending at —00 and 
encircling the origin anticlockwise, while the other reads 

^c+ioo ^ Q 


1 pC-\-lC)0 ‘ H -1 

n U"-+»> n ruT 


+ 5 )’ 


(4.3) 
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where (p{v; s) is given in (I1.18|) and c > — min{uo, ui, • ■ •, Ur}- It is worth stressing 
that '^q{u;x) can be expressed as a single contour integral 


'Sq{u;x) = 


1 1 
r(uo + 1) 27rz 


rjdw{—x) “iFi(w;uo + 1 ;m) 

q r 

X r(w) r(^/ +1 - w) 


1=1 


1=1 


r(ui +1 - w) 


, (4.4) 


where 7 encircles all non-positive integers such that Re{w} < min{/ri-|-l,..., fiq+1} 
for any w € This is a nice analogue of the definition of the function ^q{v; y) and 
can be derived as follows. First, the power series expansions for the two functions 
and qF’i (uq + 1; —xu) give us the following relation 


1 

e'"oFi{vo + 1; -xu) = X! “ ^’*^0 + 1;“); (4-5) 

fc =0 

from which, together with the definition of the function (14.21) . by term-by-term 
integration we then read off 




1 

r(uo -f 1) 


^ 2;'"'i.P'i( - + F,u) 

/c=0 


xHlfa + l + Un ( \ ) . (4.6) 

1=1 1=1 ^ ‘ > 

With this, noting that the integrand for the integral (lOl) has simple poles at 
0,-1,—2 ,..., we thereby apply the residue theorem to get the desired result. 


Proposition 4.1. Let Y be defined by and suppose that all eigenvalues 

oi,..., qm of A*A are positive. The eigenvalue PDF of Y*Y can be written as 

'Pn{xi,...,xn) = ^det[iFAr(a;i,a;j)]5=i (4.7) 

with correlation kernel 


KN{x,y) 


1 f 1 ^ 

-— / du dvu''°e~'^~^'"'^>r(u;x)^r(v,y) -TT 

Sttz Jo Jc ^ r=i 


ai 


ai 


(4.8) 


where C is a counterclockwise contour encircling —Oi,..., —ojv but not u. 


Proof. Starting with the eigenvalue PDF (12.ip of (Go -f A)*{Go -f A), application 
of [iQl Corollary 2.4] r times in succession shows that the eigenvalue PDF of Y*Y 
is proportional to 

det[r7i(a;j)]5=i det[^i(a;j)]f^.=i, (4.9) 

where rji{x) = and with T = ti ■ ■ -p 


fiix) = 


r(uo -I-1) 7(0,1)' 


dtHf^^-^l-L) 


1^1 






(4.10) 

Next, we proceed as in the proof of Proposition ll.il Our first task is to compute 
9i,j /o°° For this purpose, we note that application of the Mellin 
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transform gives 



2 /* = (j)[-aj\s)W_B{vi ui), 

Z=1 


(4.11) 


where the notation B(a,b) refers to the gamma function evaluation of the beta 
integral and s) is given in (11.181) . while use of the inverse Mellin transform 
gives 


0(y) = ^r-(-aj;?/)]^r(/ri - Pi), (4.12) 

1=1 


where is defined in ()4.3I) with q = r. Combining (j4.11[) . (jl.l7() and ()2.1ip . we 
obtain 


9^,3 = (* - l)!e“"i)'”i(-aj)]^B(ui +f,/ii - vi). (4.13) 

1=1 


According to Proposition 12.21 with G = and C = [G ^)‘, the entries 

Cij of C then satisfy 


N r 

e“'“ ^(f - l)\ L'^'l^{-ak)W_B{vi + - u/) Cjj = 5j^k- (4-14) 

i=l 1=1 


Without loss of generality we assume that ai,...,aN are pairwise distinct. The 
above equations imply 


N 




1^1 


N 

-p-r -u - ai 


(4.15) 


which can be verified by noting that both sides are polynomials of degree A^ — 1 in 
u which are equal at N different points. Using this implicit formula for {cij} and 
the integral representations 


r(z) 




dt, 


1 


f w-^e^dw, 

27ri 


0 


7 


(4.16) 
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we have from (j2.9ll that with T = ti ■ ■ - tr and W = wi ■ ■ ■ Wr 

N r « 


KN{x,y) = 


1 




(27ri)’' T { fj.i - vi ) Jr 

r r 

+ i , Yi - vi ) Cij 


dt J dw[xT/Wy 


1=1 


1=1 


N 




i=i 


(0,oo)^ 


1^1 


pxTjW 

' ' ’ r(uo + 1) 


^ pOO 

^J-Wb{vi + i,Yi - / duu''°L''°_^{u)e~'^oFi{vQ + l\-uxT/W) 

A — 1 ; 1 0 


/ = 1 
N 


(2'Ki)'' 


■^$r(-aj;y) / dt dwY[(t 


t=i 


(0,oo)^ 


^xTjW 

w, -• -e"-' -7T X 


i 


1^1 


r(uo +1) 


/•oo _ _ 

/ duu'^°e~'^oFi{vo + 1; —uxTjW) e““^ TT- 

Jo iyj " 


-u -jH 

'J ~ 


(4.17) 


Here the formulae (12.121) and (14.151) have been used in the second and third equalities 
respectively. 

Finally, recalling (14.21) we can rewrite (14.171) as 

pOO 

KN{x,y)= / (iuu''‘’e““^'r(R;a;)^$r-(-aj;y)e““^TT—(4.18) 

Jo a, ai 

If we recognise the sum over j as the sum of the residues at {a;} of the v- function 


^r{-v]y)- 


1 


N 

n 


-u - ai 


—u — V V — ai 

1=1 


by changing v to —v we then arrive at the desired result. 


(4.19) 

□ 


At this stage it would be possible to develop the theory of the corresponding av¬ 
eraged characteristic polynomials and their reciprocals, and then proceed to analyse 
their hard edge limit; recall Sections 12.21 and 13.21 However we pass on this, and 
instead analyse the hard edge phase transition analogous to the workings in Section 
13.11 Specifically, taking iV —>■ oo, we keep all Vj fixed and simultaneously let some of 
Yi,..., Yr go to oo. Without loss of generality, we suppose that for some 0 < g < r 
all ui,..., Ur; 7 * 1 ) ■ • ■; 7*g are constants, and moreover 

/iq+i,...,/ir —>■ oo as N ^ OO] (4.20) 

see [ini Theorem 2.8] for the assumptions. 

Theorem 4.2 (Subcritical kernel). With the kernel (14.81) . for a fixed nonnegative 
integer m let 


Qj = Naj, j = 1,... ,m and ak = bN, k = m + 1,..., N, 


(4.21) 
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where 0 < 6 < 1 and cti, ..., am > 0. Set cat = (1 — b)N^q+i ■ ■ ■ ^r- Under the 
assumption (I4.20|) we have 



Proof. Substituting u,v by uN,vN respectively in (14.81) . we obtain 


1 

cn 



Cn ’ 



2cAr7rz u — V 


fj (u + 6)(u + gj) 

(u + 6)(u + aj) 


{NuY'^'i'riNu-, —)^riNv; 

Cn Cn 


(4.23) 


where f{z) = z — log(6 + z). 

We can complete the proof in much the same way as in that of Theorem 13.11 
But this time we need to estimate the large N leading terms of the functions dt and 
$. That is, we have to rescale variables tq+i ,... ,tr a-nd rewrite them according to 




Cn (27ri)’’ r(uo + 1) J{o,ooy Jr 


dt dw X 


i=9+i 


q r 

exp I 

1=1 1=1 


f h-'-tq 
(1 — b)N wi - ■ -W; 


■| qFi ^uq+I; 


ju tl---tq 

1 — b wi ■ ■ ■ w.. 




(4.24) 


and 


1 rc+ioo 

^r{Nv]—) = — ds{- -—) V(7Vn;s) 




+®)n 


UifJ-i+s) 


, (4.25) 


then apply the saddle point analysis (see e.g. [60]) to the integrals over t^+i, ... ,tr 
in dtr near the saddle point to = 1, or expand the integrand in by the Stirling 
approximation formula as pLq+i ,..., —>■ oo. Tracking the same contour deforma¬ 
tions and following almost the same analysis as in Theorem 13.11 the proof will be 
done. We leave the details to the reader. □ 


The limiting kernel on the RHS of (I4.22L with the parameter uq absent and 
r -I- 1 replaced by r first appeared in [lOl Theorem 2.8] as the hard edge correlation 
kernel for a product of truncated unitary matrices. Clearly, it reduces to the Meijer 
G-kernel HU) in case q = 0. More generally, as remarked in [30] (cf. eqns (2.37) 
and (2.38) therein), it can be interpreted as a hnite rank perturbation of (11.61) . 

For the critical regime, tracking the same contour deformations and following 
almost the same analysis as in Theorem 11.31 as for the proof of Theorem 14.21 the 
required working to establish the following theorem can be given. 
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Theorem 4.3 (Deformed critical kernel). With the kernel (|4.8|1 . for a fixed non¬ 
negative integer m let 

Oj = VNctj, j = 1, ... ,m and Ofc = A^(l — t/\/N)~^, fc = m + 1,...,TV, (4.26) 

where r € R and cti, ..., a^n > 0. Set • • • Mr- Under the assumption 

(I4.20|) we have 


N^oo CN ^CN CN^ Jq J-c-ioo ^ ^ ^ 


u\^o e 

V ' 


u — V 


n 


u + a-j 


■i 


-hi, 


-hq 




- - , 'j 

3 = 1 ■' 

where 0 < c < min{(Ti,..., Om}- 




We remark that the kernels on the RHS of (14.271) reduce to the deformed critical 
kernels Kffif in (11.271) in case q = 0. These are the most general form of critical 
kernels that we have derived in the present paper. Moreover, they are new except 
for the simplest case q = r = m = 0, which as previously remarked corresponds to 
non-intersecting squared Bessel paths and has been studied in [22l EH |44] . 

As to the supercritical regime where 6 > 1, when r > 1 we have a similar 
expectation on the scaling limit (see ean. ()3.33p and relevant description below it), 
which can be stated as follows. With the kernel (14.8L for a fixed positive integer 
m let 


lim —J 
N^OO Cn 


-1), 

j = 

1,...,' 

m and ak = 

CTl,... 

7 

> 0. 

Set Cn = 

we have 



(A, 

-) 

1 

[ du [ 

\CN 

CN 3 

2'Ki 

Jo Jc 


(4.28) 


/ {-r 

Jc « 


^-U + V 


U — V 


X ^'^(m/k; K^) $q(u/K; Kjy) 


3 = 1 


u -I- a j 
v + aj 


(4.29) 


where k = {b — l)/b, ikq, are given by (14.2p . (14.31) . and C is a counterclockwise 
contour encircling —cti, ..., —cr^ but not u. 


5. Asymptotics for large parameters and variables 

5.1. Limits for large parameters. The behavior of the critical kernel (|1.24l) 
for large values of the parameters will be discussed, one of which is the confluent 
relation between correlation kernels. The first to be considered is when some of 
131,..., Ur, say Um+i,...,Ur, go to infinity. 

Proposition 5.1. Let/C*'’'’(^, ly; r) be the critical kernel (11.241) . If0<m<r, then 
as Um+i ,..., Ur ^ oo wc have 


{Um+l ■ ■ ■ Ur)]C^’^ {{Um+l ' ' ' Ur)x, (Um+1 ’ ’ ’ Ur)y,T) -S> {x, y, t) . (5.1) 

Proof. This immediately follows from the identity (13.121) for GJ ’°_|_2 and the defini¬ 
tion (11.71) for 
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The above confluent relation allows for a natural interpretation, particularly 
in the original finite matrix dimension. Actually, in ()1.4|) substituting all Gj as 
square matrices being distributed according to the joint density proportional to 
det'^^ {G*Gj) exp{—tr(G'*G'j)} (see the relevant description below (11.41) 1. then a sad¬ 
dle point approximation shows that all Gm+i, ■ • •, Gr go to the identity matrix of 
order N as Um,. ■ •, Ur oo. Thus these matrices do not contribute to the hard 
edge state. 

A similar effect happens in relation to the parameter uq associated with Gq, 
although now we find that a different rescaling is necessary, and furthermore that 
the limiting kernel is now subcritical. 


Proposition 5.2. Letry; r) be the critical kernel (11.241) . Forr > 1, we have 
lim (5.2) 

Vq—^OO 

where is given by mBi)- 


Proof. Substituting u,v by ^/vqu and respectively in (11.241) . we get 

t) = ^ [ du [ dv 
2Tn Jo 


p-VQU(n)-f{v)) _ ^_ 

_ _ ^-ry/v^u-\-Ty/v^v 


U — V 


X uoG,^’° 


0, t -+2 


0 , 


i/oua;')Gd++2( „ vovy') (5.3) 

/ \ z/q , iy\ ,... , Uj’ , U / 


where f{z) = — log 2 12. 

Choose one saddle point zq = 1 from f'{z) = 0 and deform iM. as the union of 
one closed clockwise contour C encircling the interval [0,1) and the vertical line 
X = 1. Note that as uq —t oo 


IZqG, 


0x+2{ 


0, -Uo, -Pi, ■ • ■ , -Pr 


UoUX 


[ ^r+1,0 / 


^1,0 1 
^0,r+l I 


0 ,-Pi, 


PO, Pi,• ■ •, 


izovy] 


Pi,.. 


vy), (5.4) 


proceeding as in the proof of Theorem 13.11 we can show that the dominant contri¬ 
bution comes from the range of u £ [0,1) and v € C. Finally, application of the 
residue theorem gives the proof. □ 


Similarly, for the large negative r, we observe a transition from the critical kernel 
to the Meijer G-kernel. This is to be expected, as then the parameter b in (11.111) 
enters the subcritical regime 6 < 1, since effectively 6 = (1 — t/N)~^. 

Proposition 5.3. Let {f,r];T) be the critical kernel ()1.24l) . Then we have 


lim x / T ,- y / T -, T ) = K ^’''{ x , y ). 


Proof. Substituting u,v by —tu and —tv respectively in (ll.24p . we get 

, 1^0 g-PifG)-fG)) 


- x / T ,- y / T -, T ) = -^ [ du [ dv (-) 
2 tti Jo JiR \v/ 


U — V 


xG. 


1.0 / 

0,r-|-2 


g: 


r-l-l.o f 
0,r+‘2 


uo,ui,...,ur,0 


vy 


,0,-z/Q,-z/i,...,-Ur 

where f{z) = —z + z^ 12. Proceeding as in the proof of Proposition [521 the 
result follows. 


(5.5) 


(5.6) 

sought 

□ 























30 


PETER J. FORRESTER AND DANG-ZHENG LIU 


Lastly, as to the critical kernel on the RHS of (I4.27|) with m = 0, the functions 
and defined in (13.3711 and (I3.38p . there exists similar asymptotic 

behavior for large parameters as in the above three propositions, but we refrain 
from writing them down. 


5.2. Conjectures and open problems. In the concluding section of [53] a num¬ 
ber of questions, mostly relating to asymptotics, were posed in relation to the 
kernel (HSj). As we will indicate, these all carry over to the critical kernel (|1.24|) . 
It is also the case that the conjectured behaviours are all closely related to analo¬ 
gous expected asymptotic properties of the finite N kernel (|1.22|) . Two classes of 
asymptotic problems stand out. 

The first is to establish the global scaling limit of the critical one-point function. 
For this we expect 

\im N'^KNiN^+^x^N^+^x) ^-lmG{x-iO), (5.7) 

A-^oo ai=N IT 

where w{z) := zG{z), satisfies the algebraic equation 

^r-S3/2 _ ^^ 1/2 = 0. (5.8) 


The latter is known to specify the Raney distribution with parameters (3 -I- 2r, 2), 
which according to free probability theory is the global density for the matrix dEl 
in the critical case (see e.g. |3TJ Remark 3.4]). In the case of the global limit ()5.7I1 
with a; = 0 (^ = 1, ..., N), a recent achievement [48] has been the use of the double 
contour integral formula (12.211) to deduce that (15.71) with w{z) := zG{z) satishes 
the algebraic equation 

— zw + z = 0. (5-9) 


The latter specifies the Raney distribution with parameters (r -|- 2,1), also known 
as the Fuss-Catalan distribution with parameter r -|- 1 [56], and should give the 
asymptotic behavior of global density for small argument throughout the subcritical 
regime. In the supercritical regime, from a macroscopic viewpoint the number of 
random matrices in the product is effectively r, since A dominates Gq and 
moreover A is proportional to the identity. This implies that the corresponding 
asymptotic behaviour of the global density near the origin now corresponds to that 
of the Fuss-Catalan distribution with parameter r. 

To see the relevance of (|5.7I) to the asymptotics of the density in the critical hard 
edge scaled state, KA’^{x,x), we recall (cf. [ST] Cor. 2.5]) that it can be deduced 
from (15.9p that for small x the global density has its leading asymptotics given by 
(I1.12P . In keeping with the discussion in the concluding section of [29] , this should 
be the leading large x asymptotic form of ’''{x, x). Combining this with the small 
X asymptotic form (11.51) for the Fuss-Catalan density as applies to the subcritical 
and supercritical regimes (the latter with r i—>■ r — 1 as already commented), we 
therefore expect 




1 

sin 

TT 

- X 

l+7T2^ 

TT 


r + 2 


1 

sin 

2tt 

- X 

^+7+37: 

TT 


2r + 3 


1 

sin 

TT 

-7 ^ 

1+ETT^ 

TT 


r -b 1 



IC'^’^’^*\x,x) 


x—^oo 


(*) = subcritical 
(*) = critical 
(*) = supercritical. 


(5.10) 
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In general if the global density at the hard edge diverges as x~p , then the ex¬ 
pected number of eigenvalues in the interval (0, s) is proportional to N. For 
this to be of order unity we must scale s i—>■ Taking into consideration 

the scaling x i-A- N'~^^x already present in (15.71) . this suggests that the appropriate 
hard edge scalings are 

{ {l/N)x, subcritical 

{1/^/N)x, critical (5-11) 

no change, supercritical, 

in agreement with those used in the main body of the text. 

The second class of asymptotic of the type identified in the concluding section 
of [29] is to compute the leading asymptotic form of the off diagonal analogue of 
the LHS of (|5.7L namely 

kN{x,y)-.= N'^+^KN{N'^+^x,N''+^y), x^y, (5.12) 

where the dot above := indicates that terms which oscillate and average to zero 
are to be ignored. To see the interest in this quantity, note from (12.101) that the 
truncated (or connected) two-point correlation ^^(aii,X 2 ) := P(2),Ar(3ii, 2 : 2 ) — 
P{i),n{xi)p{i),n{x2) is given hy pJ^-^ j^{xi,X2) = -Kn{xi,X2)Kn{x2,xi), so knowl¬ 
edge of the asymptotics of K{x,y) tells us the asymptotics of 

pf 2 )M^, y) ■■= (5-i3) 

With G = yi^j) denoting a linear statistic in the bulk scaled system, in view 

of the formula (see e.g. [23 eqn. (14.38)]) 


VarG = dxi dx 2 g{xi)g{x 2 )pf 2 ) 

Jo Jo 

poo 

-I-/ g(^x)p(^i){N^~^^x) dx (5.14) 

JQ 

one sees that (15.131) (sometimes referred to as a wide correlator; see e.g. [36]) es¬ 
sentially determines the large N form of this fluctuation, which is expected to be 
G(l) (see e.g. [23 §14.3]). 

As a concrete example of this second type of asymptotics, consider the simplest 
case of (O, namely r = 0 and A = 0. The squared singular values correspond 
to the eigenvalues of GqGq, where Go is a {N + vq) x N standard complex Gauss¬ 
ian matrix. This class of random matrices is referred to as the complex Wishart 
ensemble (see e.g. [23 §3.2]). For this ensemble it is a known result that [TO] 


2VVf2),Ar(A^a:, IV?/) ~ 


1 1 {L/2){x + y) - xy 

27r2 (x - ?/)2 (x(L - x)y{L - y))^/^ ’ 


x^y, (5.15) 


with L = 4, and where the dot above the asymptotic sign denotes a restriction to 
non-oscillatory terms. 

Suppose now that in the definition (15.121) of Kn{x, y) we introduce a scale factor 
L and compute instead the asymptotic form of /L)K]s[{N‘^~^^x/L, N'^^'^y/L). 

For the complex Wishart ensemble the RHS of (15.151) with L a variable results. For 
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general r, if the original leading asymptotic form of y) was R{x, y), this will 

now equal (\jL?'')R{xjL^yjV). Following [TU] we expect that 


where R^{x,y) is the leading non-oscillatory large x, large y asymptotic form of 
the hard edge scaling of j^{x,y). In the context of the present setting this 
corresponds to seeking the large x, large y form of K}'’^{x,y). In the case of the 
complex Wishart ensemble, (15.161) applied to (15.151) predicts that 


h,T/ N 

P(2) i^^y) 


1 {xlyY^'^ + jylxY^^ 

{x — yY 


(5.17) 


which is in fact a known exact result (see e.g. [m eqn. (7.75)]). The analogue of 
(I5.17|) is known for the case r = 1, A = 0 of ()1.4I) [29l eqn. (5.28)], but the analogue 
of (|5.15|) is yet to be obtained. As discussed in [29], knowledge of an asymptotic 
form such as (|5.17l) is of interest for the computation of the variance of a scaled 
linear statistic at the hard edge, Ga = when a —>• oo, which is given 

by 


lim Var Ga '■= lim 

cn—>-oo a^oo 


dX\ I dX2 
Jo 


xg{Xi/a)g{X2/a)p^^^'''{Xi,X2) + J dXg{X/a)p'^{^{X)y (5.18) 

A number of challenges for future research present themselves from the above 
discussion. We conclude this section with a list of a few more. 


• Under the assumption of ai = • • • = oat = bN with b > 0, verify the sine- 
kernel in the bulk and Airy-kernel at the soft edge for (|1.22|) and (|4.8I) (see 
recent monographs [6l|2ll|27l|58] for the sine and Airy kernels and [48] for 
recent progress on the random matrix products). 

• Under the assumption of Om+i = ■ ■ ■ = aw = bN with 6 > 0, by tuning the 
parameters ai,..., am verify the BBP transition for (11.221) and (14.81) (cf. 

mm)- 

• Verify the transitions from the critical kernels (|1.24|) and (|4.27l) to the sine- 
kernel and to the Airy-kernel (cf. [27l Exercise 7.2] and [29]). 
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